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' Abstract 

The symmetry energy effects on the location of the inner edge of neutron star crusts are studied. Three 
04 . phenomenological models are employed in order to check the accuracy of the well known parabolic approximation 

of the equation of state for asymmetric nuclear matter in the determination of the transition density nt and 
transition pressure Pt . The results corroborate the statement that the error due to the assumption that a priori 
the equation of state is parabolic may introduce a large error in the determination of related properties of a 
neutron star as the crustal fraction of the moment of inertia. 

o 
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1 Introduction 

' Neutron stars (NS) are extraordinary astronomical laboratories for the physics of dense neutron-rich nuclear matter 
• [IlEj- The main parts of a NS are the crust and the core. The latter, divided into the outer core and the inner one, 
CNJ ' has a radius of approximately 10 km and contains most of the star's mass while the crust, with a thickness of about 
, 1 km and containing only a few percent of the total mass, can also be divided into an outer and an inner part. A 
■ very important ingredient in the study of the structure and various properties of neutron stars is the equation of 
. . , state (EOS) of neutron-rich nuclear matter [3,. 

One of the most important predictions of a given EOS is the location of the inner edge of a neutron star crust. 
The inner crust comprises the outer region from the density at which neutrons drip-out of nuclei, to the inner edge 
5J] ' separating the solid crust from the homogeneous liquid core. At the inner edge, in fact, a phase transition occurs 
from the high-density homogeneous matter to the inhomogeneous one at lower densities. The transition density 
takes its critical value Uc when the uniform neutron-proton-electron matter (npe) becomes unstable with respect 
to the separation into two coexisting phases (one corresponding to nuclei, the other one to a nucleonic sea) [3]. 

While the density at which neutrons drip-out of nuclei is rather well determined, the transition density nt at the 
inner edge is much less certain due to our insufficient knowledge of the EOS of neutron-rich nuclear matter. The 
value of nt determines the structure of the inner part of the crust. If sufficiently high, it is possible for non-spherical 
phases, with rod- or plate-like nuclei, to occur before the nuclei dissolve. If nt is relatively low, then the matter 
undergoes a direct transition from spherical nuclei to uniform nucleonic fluid. The extent to which non-spherical 
phases occur will have important consequences for other properties determined by the solid crust [4 . 

In general, the determination of the transition density nt itself is a very complicated problem because the 
inner crust may have a very complicated structure. A well established approach is to find the density at which 
the uniform liquid first becomes unstable against small-amplitude density fluctuations, indicating the formation of 
nuclear clusters. This approach includes the dynamical method 'SJ [SI [71 |H1 [HI UHl IHl [H] , the thermodynamical 
one [31 [131 [m [H] , and the random phase approximation (RPA) [16. [T7] . 

Theoretical studies have shown that the core-crust transition density and pressure are very sensitive to the 
density dependence of the nuclear matter symmetry energy [51[TTlP^[T5 1 [n i [T51[TO l [r7 l [T51ptl[ ^ l ^ |^[ ^ . 
However at present, the symmetry energy is almost unknown experimentally in the density regime appropriate for 
the interior of neutron stars |24[ I25[ 126] . In spite of the experimental uncertainty of the symmetry energy, there 
are many theoretical considerations of the symmetry energy categorized mainly in phenomenological and effective 
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field theoretical approaches. The aim of this work is: a) to apply a momentum dependent interaction model (MDI) 
as well as two additional non relativistic models (Thomas-Fermi and Skyrme type) to determine the transition 
density and transition pressure corresponding to the edge of a neutron star crust and b) to check the accuracy of 
the parabolic approximation widely used in the literature and applied to neutron star research. 

The article is organized as follows. In Sec. 2 we review the Taylor expansion of the energy while in Sec. 3 we 
present the thermodynamical method for the determination of the transition density and pressure of the inner edge 
of a neutron star crust. The nuclear models employed in the the present work are presented in Sec. 4. The results 
are presented and discussed in Sec. 5. Finally Sec. 6 summarizes the present study. 



2 Taylor expansion of the energy 

The energy per particle E{n^ I) in cold asymmetric nuclear matter can be expanded around / = as follows 

E{n, I) = E{n, / = 0) + EsymAn)!^ + Esy^A{n)I^ + ... + E.ymMn)!^'" + ■■■, 



(1) 



where the total baryon density n = jip + Un, Up (rip) is the proton (neutron) density, / = (n„ — np)/n is the 
asymmetry parameter and E{n, / = 0) is the energy per baryon of the symmetric nuclear matter, while the 
coefficients of the expansion are 



Esym,2{n') = 
Esym.2k{n) 



1 d^E{n,I) 



2! dP 



E. 



symA 



(n) = 



1 d^Ein.I) 



1=0 



4! 9/4 



7=0 



1 d^''E{nJ) 



{2k)\ dP'' 



(2) 



1=0 



In II]) , only even powers of / appear due to the fact that the strong interaction must be symmetric under exchange 
of neutrons with protons i.e. the contribution to the energy must be independent of the sign of the difference 
nn — np. The second order approximation of the expansion ([Ij is written as 



E{n, I)^E{n,I = 0)+ Esyr,-,,2 (n)!^ 



(3) 



We expect, from a mathematical point of view, that the above expansion is accurate at least close to 7 = (the 
case of symmetric nuclear matter). However for the majority of the energy functionals, the above approximation 
works well for higher values of the asymmetry parameter / and even close to the value 1=1 corresponding to 
the pure neutron matter. Thus, in cases for which the expansion, in a good approximation, is independent of the 
asymmetry parameter /, the symmetry energy can be defined as 



Esy„,{n) = E{n, 1 = 1)- E{n, I = 0), 

and the energy per baryon is written 

E{n, I) = E{n, I = 0) + {E{n, 1 = 1)- E{n, I = 0)) P. 



(4) 



(5) 



The knowledge of the equation of state of neutron rich matter is fundamental in astrophysical applications. For 
example it is the basic ingredient for the study of /3-stable matter characteristic for the interior of neutron stars. 
Actually, for the most of the equations of state (mainly those originating from microscopic calculations or those 
coming from relativistic mean field theories) only the energy of symmetric nuclear matter and pure neutron matter 
are determined and the definition of the symmetry energy from equation Q is almost unavoidable. The question 
naturally arising is the magnitude of the width of the error introduced by assuming a priori that the equation of 
state is parabolic according to relations ([T]) and ([5]). 

In the present work we define in addition the second order expansion of the form 



E{n, I) ~ E{n, I = 0) + E,y„,^2 (n)P =E{n,I = 0) + 



1 d^E{n,I) 



2! dP 



(6) 



which is similar to expansion ([5]) replacing the quantity Esym{ri) = E{n,I = 1) — E{n,I = 0) by Egymfiin) = 
ld^E{nJ) 



1=0 



2! dP 

Before trying to check the accuracy of the approximation ([5]) it is worth to compare the density dependence of 
the above two definitions of the symmetry energy i.e. Esym.ai''^) Esy,n(n) given by ([2]) and Q respectively. 
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3 The thermodynamical method 



The core-crust interface corresponds to the phase transition between nuclei and uniform nuclear matter. The 
uniform matter is nearly pure neutron matter, with a proton fraction of just a few percent determined by the 
condition of beta equilibrium. Weak interactions conserve both baryon number and charge [3] , and from the first 
law of thermodynamics, at temperature T = we have 

du = —Pdv — /idq, (7) 

where u is the internal energy per baryon, P is the total pressure, v is the volume per baryon { v — 1/n where n is 
the baryon density) and q is the charge fraction (q = x — Ye where x and Y,, arc the proton and electron fractions 
in baryonic matter respectively). In /3-equilibrium the chemical potential fi is given hy ft — fj,n — fJ-p = fJ-e where 
/Ltp, Hn and fie are the chemical potentials of the protons, neutrons and electrons respectively. The stability of the 
uniform phase requires that u{v,q) is a convex function [27 . This condition leads to the following two constraints 
for the pressure and the chemical potential 



dp\ _ zap 

dv J „ \ dq ) „, \dv 



> 0, (8) 



dfi 



> 0. (9) 



It is assumed that the total internal energy per baryon u{v, q) can be decomposed into baryon (i?jv) and electron 
[Ee) contributions 

u{v,q) = EN{v,q)+Ee{v,q). (10) 

The related theory has been extensively presented in our recent publication [20.^. We consider the condition of 
charge neutrality q = Q which requires that x = Ye- This is the case we will consider also in the present study. 
Hence, according to Ref. [3D] the constraints ([SJ and (0), in the case of the full EOS (FEOS) and the parabolic 
approximation (PA), are written 



^ dE{n,x) ^d^E{n,x) fd^E{n,x) Yfd^E{n,x)\ ' ^ 



CpA(n) 



^d'^E{n,x = 0.5) 

dn^ dn 
^ d'^Esymjn) , ^ dEsymin) ^ 



dE(n,x = 0.5) , 
In + (1 - 2x)^ 



(12) 



1 



dv? 



dn 



Esyra (^) 



^dE gyyyi{Tl f 

dn 



> 0. 



For a given equation of state, the quantity Cfeos(?^) (or Cpa('^)) is plotted as a function of the baryonic density 
n and the equation Cfeos(") — (or C-px{n) = 0) defines the transition density n^. However, what remains is the 
determination of the proton fraction x (as a function of the baryon density n) in /3-stable matter. In this case we 
have the processes 

n — > p + + Pe p + — > n + Ve (13) 

which take place simultaneously. We assume that neutrinos generated in these reactions have left the system. This 
implies that 

A = A'n - = A^e, (14) 

Given the total energy density of the baryons eat = e(ri„,rip), the neutron and proton chemical potentials can be 
defined as (see also Ref. [28 ) 



den 
dnn 



dn 



N 



It is easy to show that after some algebra we get 



d{eN/n) 
dx 



dEN \ 
' dx J, 



(15) 



(16) 



The charge condition implies that Ue = np = nx or kp^ = kp (where kp are the fermi momenta). In addition, the 
chemical potential of the electron is given by the relation (relativistic electrons) 



Me 



k'^p c? + TO^c* ~ kp^c — hc{3n^nx) 



1/3 



(17) 
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Finally, from equations ([TB]) and (fTT]) one has 



dE 



N 



dx 



(18) 



Equation ([18]) is the most general relation that determines the proton fraction of /3-stable matter. In the case of 
the PA the above equation, with the help of Eq. ^ is written as 



4(1 - 2x)Esy,n{n) = hci'iTi'^nef/^ = hc{'iTi'^nxY/^. 



(19) 



The pressure Pt at the inner edge is an important quantity directly related to the crustal fraction of the moment of 
inertia, which can be measured indirectly from observations of pulsars glitches 0] . The total pressure is decomposed 
also into baryon and lepton contributions 



where 



P{n,x) = PN{n,x) + Pe{n,x), 
,dEN 



PN(n, x) = n 



dr 



The electrons are considered as a non-interacting Fermi gas. Their contribution to the total pressure reads 



Pe{n,x) = 



1 



he 



(Sn'^xrij 



4/3 



127r2 (?ic)3 127r2 

The transition pressure, in the case of the FEOS, is given now by the equation 



P 



FEOS 



(nt,xt) = n\ 



dE 



N 



dr 



he 
12^ 



(Stt xtut 



,4/3 



(20) 
(21) 

(22) 
(23) 



where Xt is the proton fraction related to the transition density. In the case of the PA Pf is given by the relation 



2 / dE{n,x^ 0.5) 



\ dn 

/ 2 \4/3 

[Stt XtUt) 



dEsymin) 



{l-2xt)' 



(24) 



3.1 An application: crustal fraction of the moment of inertia 

The crustal fraction of the moment of inertia AI/I can be expressed as a function of M (star's total mass) and R 
(star's radius) with the only dependence on the equation of state arising from the values of Pt and rit. Actually, the 
major dependence is on the value of Pt, since rit enters only as a correction according to the following approximate 
formula [29] 

A/ _ 287r Pti?3 (1 - 1.67/3 - 0.6/3^) / 2Pt (1 + 7P){1 - 2P) 
~T ~ 3Mc^ ^ V ^ ntmc^ /P 

where /? — GM/Rc^. The crustal fraction of the moment of inertia is particularly interesting as it can be inferred 
from observations of pulsar glitches, the occasional disruptions of the otherwise extremely regular pulsations from 
magnetized, rotating neutron stars [Tl]. Link et al. [^ showed that glitches represent a self-regulating instability 
for which the star prepares over a waiting time. The angular momentum requirements of glitches in the Vela pulsar 
indicate that more than 0.014 of the star's moment of inertia drives these events. So, if glitches originate in the 
liquid of the inner crust, this means that AI/I > 0.014 




4 The models 

In the present work we employ three different phenomenological models for the energy per baryon of the asymmetric 
nuclear matter having the advantage of an analytical form. The MDI model has been extensively applied for neutron 
star studies, can reproduce the results of more microscopic calculations of dense matter at zero temperature and 
can be extended to finite temperature [301 ISIl 1311 1331 1311 133 ■ The Skyrme model using various parametrizations 
can also be applied both in nuclear matter and in finite nuclei [36[ 137) . Finally we employ a version of the Thomas- 
Fermi model, which was introduced by Myers et. al. |38| and has also been applied for the study of finite nuclei 
and in a few cases for high density nuclear matter applications [39) . 
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4.1 MDI model 



The model used here, which has aheady been presented and analyzed in previous papers [SUl 1311 1321 [Ml IMl ISU : is 
designed to reproduce the results of the microscopic calculations of both nuclear and neutron-rich matter at zero 
temperature and can be extended to finite temperature |30) . The energy density of the asymmetric nuclear matter 
(ANM), in MDI model is given by the relation 



e(n„,np,r = 0) = efc,„(n„,T = 0) + el^Jnp,T = 0) + K:n^(?^T^, "p, T = 0), 



(26) 



where the first two terms are the kinetic energy contributions on the total energy density, while the third one is 
the potential energy contribution. The energy per baryon at T = 0, is given by 



Ein,I) = ^E°pU^^^[il + lf/^ + {l-lf/^ 
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1/3^ 



(27) 



In Eq. (j27p . I is the asymmetry parameter (/ = (ri„ — np)/n) and u = n/riQ, with uq denoting the equilibrium 
symmetric nuclear matter density, uq — 0.16 fm~^. The parameters A, B, cr, Ci, C2 and B' which appear in the 



description of symmetric nuclear matter are determined in order that E{n - 
fm"'^, and the incompressibility is K — 240 MeV and have the values A 
Ci — —83.84, C2 = 23 and B' = 0.3. The finite range parameters are Ai — 
Fermi momentum at the saturation point rig. 



no) — mc — 
= -46.65, B 
l.bkp and A2 



16 MeV, no = 0.16 
: 39.45, a = 1.663, 
: 3fc^ and kp is the 



Table 1: The parameters for neutron rich matter of MDI model. 



F{u) 


xo ^3 Zx Z-i 


\fu 


0.376 0.246 -12.23 -2.98 


u 


0.927 -0.227 -11.51 8.38 


2wV(l + m) 


1.654 -1.112 3.81 13.16 



The additional parameters xq, a^s, Z\^ and Z2 employed to determine the properties of asymmetric nuclear 
matter are treated as parameters constrained by empirical knowledge |30| and presented in Table 1. By suitably 
choosing the parameters a;o, xa, Zi, and Z2, it is possible to obtain different forms for the density dependence of 
the symmetry energy Esymiu}- The nuclear symmetry energy is parametrized according to the following formula 

Esyrr.{n) = 13^^/3 \lF{u), (28) 

where the first term of the right-hand side of Eq. ((28)) represents the contribution of the kinetic energy and the 
second term is the contribution of the interaction energy. In general, in order to obtain different forms for the 
density dependence of Esym{n), the function F{u) can be parameterized as follows |30] 

F{u) = ^, F{u)^u, F{u) ^2u^ /{l + u). (29) 

In the present work its parametrization corresponds to the models MDI-1, MDI-2 and MDI-3 respectively. Numer- 
ical values of the parameters that generate these functional forms are given in Table 1. It is worthwhile to point 
out that the above parametrization of the interaction part of the nuclear symmetry energy is extensively used for 
the study of neutron star properties j30l [28] as well as the study of the collisions of neutron-rich heavy ions at 
intermediate energies [40l [41] . 

The pressure of the baryons, at T = 0, defined as 

on 
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with the help of Eq. (|27p takes the analytical form 



Pin, I) = ino£;Ou5/3 [(1 + 7)5/3 _^(^_j)5/3- 



(31) 
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4.2 Skyrme model 



The Skyrme functional providing the energy per baryon of asymmetric nuclear matter is given by the formula 

MM 



3 ?iV / 3^2x2/3 



10 m 



^ ) n^^^F^/sil) + hnn [2{xo + 2) - (2x0 + 1)^2(7)] 



+ j^hn-+' [2{x3 + 2) - (2x3 + 1)^^2(/)] 



(32) 



3 /3^2 



40 V 2 
1 



2/3 



.5/3 



(ti(xi+2) + <2(a;2 + 2))f^5/3(J) 



+ -ihi2x2 + l)-hi2xi + l))Fs/3iI) 



where F™(/) = - [(1 + /)™ + (1 - /)"] and the parametrization is given in Refs [551157] . 
4.3 Thomas-Fermi model 

The model's prediction for the equation of state reads E{n,I) = Tor]{u,I), where [55] 

r]{u,I) =a{I)n^ -b{I)n^ +c{I)n^, ^2=(n/no)^^^ no = 0.16114 fm"^ 



(33) 



riQ and Tq are the saturation density and Fermi energy of standard nuclear matter accordingly as predicted by the 
model and the coefficients a{I), b{I), c{I) are the following functions of / 



2(i-7i)(p» + <;°)-T. 



{5p q — q ) for n„ > Up 
{5p^q^ — p^) for n„ < iip 



(34) 



&(/) = - [a/ (/ + g^) + 2a„pV], 



c{I) = ^[Bi{p^ + q^) + B^p^q^p^ + q^)] 



(35) 
(36) 



where p = (1 + lY^^, g = (1 — 7)^/'^. The interaction strengths on, Uu, Bi, have the following values: 

7j = 0.25198, 7„ = 0.88474, ai = 0.7011, a„ = 1.24574, = 0.22791, S„ = 0.8002. 



5 Results and Discussion 

In Fig. 1 we compare the density dependence of the quantities Esym,2{n) 



1 d^E{n,I) 



2! dP 



and Esy„i{n) 



1=0 



E{n, I — 1) — E{n, I ~ 0). Both definitions of the symmetry energy exhibit a similar behavior for low values of the 
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Figure 1: The symmetry energy Esym{u) defined by Eq. (jH) versus the symmetry energy Esym,2{u) defined by 
Eq. ([2|) as a function of the density fraction u = n/riQ for the three considered models. 




Figure 2: The density dependence of the ratio EgymA^'^) I ^sym.2{u) for the three considered models. 




Figure 3: The density dependence of the proton fraction x obtained from the considered models employing the 
FEOS well as the PA. 
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Figure 4: The derivative of the symmetry energy dEsym,2{n) / dn versus the transition density nt corresponding to 
the considered models a) employing the full EOS and b) employing the parabolic approximation. 



baryon density (up to u ~ 1.5). The difference, which is model dependent, is more pronounced for higher values of 
n. This is a direct indication that the two approaches affect the equation of state of high density nuclear matter at 
the inner core of neutron stars. It is worth noticing that for the three MDI models and the three Skyrme models 
the symmetry energy is an increasing function of the density. In TF model the symmetry energy is increasing an 
function of u only up to m ~ 1.5 and then decreases rapidly with u. This has a dramatic effect on neutron stars 
applications mainly on those connected with the high density equation of state. 

In order to check the accuracy of the parabolic approximation we display in Fig. 2 the density dependence of 
the ratio Esym,i{u) / Esym,2{u). It is seen that in most cases the contribution of the fourth-order term Esym,iiu) is 
less that 4% compared to the second order one Esym,2{u)i at least for low values of the density {u < 2). However, 
in four of the considered models, i.e. the MDI-1, MDI-3, the SKa and TF the contribution of Esym,A:{u) increases 
rapidly with the density and consequently influences the high density equation of state and this must be taken into 
account. In the specific case of the TF model the contribution of Esym.iiu) becomes comparable to Esym,2{u) even 
for low values of u, making the parabolic approximation problematic. 

In Fig. 3 we compare the density dependence of proton fraction x in /3-stable nuclear matter determined by 
employing the full EOS and the parabolic approximation. In the case of the MDI and Skyrme models there is a 
difference (depended on the specific model) mainly between 4 % and 20% for low values of the baryon density (up 
to M ~ 0.7). However, the parabolic approximation is good for higher values of n. We expect that the above density 
dependence of x will be reflected also on the values of the transition density nt and pressure Pt- In the case of 
the TF model only for low values of the density {u < 1) the two approximation produce similar results. Moreover, 
it is concluded that the parabolic approximation does not affect appreciably the onset of the URCA process with 
critical values xurca — 11% (without muons) and xjjrca — 14% by including muons. 

In Fig. 4(a), we present the correlation between the derivative of the symmetry energy Esym,2{ii') (for n = l/2no, 
n = Uq and n = 3/2no) for the MDI and Skyrme models and the transition density nt in the case of the FEOS. 
The most distinctive feature is the concentration of the data from the case n ~ l/2no (where they are rather 
random) to an almost linear relation in the case n = S/2nQ. To further illustrate this point, we plot in Fig. 5(a) 
the derivative of the symmetry energy E'^^^ 2(^) at the baryon density n — 3no/2 versus nt- A linear relation is 
found of the form 

KymA^na/2) = 400.18 - 3483.3^*. (37) 
In addition in Fig. 5(b) we plot nt versus E'^y^-^ jC^'t-o/S) with the linear correlation 

nt = 0.11372 - 0.0002785i;^j,„^2(3no/2). (38) 

We find that the values of nt depend on the trend of the equation of state and the symmetry energy depends not 
only on values close to the saturation density no but even on higher ones. It is of interest to see that, according to 
Fig. 4(b) in the case of the PA the data are scattered even for values close to n = 3/2no. It is presumed that the 
use of the PA may affect the linear correlation between i?^j,„j 2(") and nt for n — 3/2rio. 

In Table 2 we present the values of nt and Pt determined by employing the full EOS and the parabolic approx- 
imation for the models considered in the present work. In most of the cases (the only exception is the Sly4 case) 
the use of the PA increases the values of nt by 10-15 % or even more (see the case MDI-3). The effect of the PA is 
even more dramatic in the case of Pt . In most of the cases Pt increases significantly compared to the FA (the only 
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n,(fm") dE^^^^(n)/dn(n=3/2nj 

Figure 5: (a) The derivative of the symmetry energy dEsym,2{n — 3/2nQ)/dn versus the transition density, for 
the FEOS, corresponding to the considered models and (b) rit versus dEsym,2{n = 3/2no)/dn. The sohd hues 
correspond to the least-squares fit (FIT) expressions ([57]) and psp respectively. 

Table 2: The transition density nt (in fm~^) and pressure Pt (in MeV) obtained from the considered models by 
employing the full EOS as well as the parabolic approximation. 



approach 


MDI-1 MDI-2 MDI-3 TF Sly4 SKI4 Ska 


nt (FEOS) 


0.0771 0.0706 0.0508 0.1368 0.0979 0.0806 0.0788 


nt (PA) 


0.0866 0.0934 0.1027 0.1419 0.0942 0.0909 0.0933 


Pt(FEOS) 


0.3456 0.2338 0.3133 2.9020 0.5781 0.3365 0.5295 


Pt(PA) 


0.2975 0.7055 1.2482 3.7300 0.5459 0.4965 0.8651 



exceptions are the MDI-1 and the Sly4 models). The increase is even two or three times. In order to clarify further 
this point we plot in Fig. 6 the values of Pt versus nt for the considered models. We see the strong dependence of 
Pt on ut in the case of PA compared to FEOS. The above results indicate that one may introduce a large error by 
employing the parabolic approximation in order to determine the value of Pt- 

The values of Pt, according to Eq. ((25|) . strongly influence the value of the crustal fraction of the moment of 
inertia. Actually, the crustal fraction A/// depends not only on Pt and nt but on the mass M and radius of R of 
the neutron star as well. Actually M and R depend also on Pt and nt . In order to have a qualitative picture of the 
effect of Pt and nt on AI / 1 we can consider the maximum mass and the corresponding radius R of the neutron 
star coming from the MDI modelI12] and neglecting the effects of Pt and nt (which is a good approximation). We 
determine the ratio A/// by considering the full EOS and the parabolic approximation. The results are presented 
in Table 3. Its obvious that there is a strong difference between of two approaches. For example in the case of 
MDI-2 the PA increases more than two times the value of AI/I compared to the FEOS, while in the case of MDI-3 
the increase is even higher. 

Table 3: The crustal fraction of the moment of inertia for the three MDI models employing the full as well as the 
parabolic approximation. 



Model 


M,nax 1 Mq 


R{Km) 


A///(FA) 


A///(PA) 


MDI-1 


1.895 


10.112 


0.0076 


0.0068 


MDI-2 


1.935 


10.570 


0.0062 


0.0163 


MDI-3 


1.952 


10.933 


0.0086 


0.0287 



Finally, the limit AI/I = 0.014 constrained by Link et.al.I^nj limits also the masses and the radii of the neutron 
star for specific values of Pt and nt. It is concluded that the strong dependence of AI/I = 0.014 on Pt puts also 
strong constraints on the allowed pairs of M and R. In Fig. 7 we plot the M-R constrained relation for the Vela 
pulsar where AI/I > 0.014 obtained from the full and parabolic approximation in the case of the MDI model and 
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Figure 6: The transition pressure Pt versus the transition density rit employing the fuh and the parabohc 
approximation. 




Figure 7: The M-R constrained relation for the Vela pulsar where AI/I > 0.014 obtained from the FEOS and PA 
in the case of a) the MDI model and b) the Skyrme model. The constraint implies that allowed masses (in M0) 
and radii lie to the right of the line. The two additional constraints are taken from Ref. [H] {rit — 0.065 fm^'^ and 
Pt = 0.26MeV fm-3) and Ref. [H] (m = 0.075 fm^^ ^^^^ ^ 0.65MeV fm~3). 



the Skyrme model. It is obvious that the implication of the full EOS imposes more restrictive constraints compared 
to the parabolic one. This effect is more pronounced in the case of the MDI models compared to Skyrme model. 
Consequently, the use of the PA may also introduce a large error in the determination of the minimum radius (for 
a fixed value of a mass) of a neutron star. 



6 Summary 

The study of symmetry energy effects on location of the inner edge of neutron star crusts is the main topic of the 
present work. We employ three different phenomenological models for the prediction of the EOS and we determine 
in the framework of thermodynamical method both the transition density and transition pressure corresponding to 
the inner edge of a neutron star crust. Actually, rit and Pt have been determined separately by employing the full 
equation of state of asymmetric matter and its parabolic approximation. Our results confirm the statement that 
employing the PA instead of the full EOS may influence appreciably the values of rit and even more the values of Pt. 
We found also that the implication of the full EOS imposes more restrictive constraints, on the M-R constrained 
relation for the Vela pulsar, compared to the parabolic one. Consequently, the error introduced by assuming that 
EOS is parabolic may induce also a large error in the determination of related properties of a neutron star as the 
crustal fraction on the moment of inertia. Finally, we found that there is a linear correlation between the derivative 
of the symmetry energy E'^y^ jl*^) at the baryon density n — 3no/2 and the transition density nt for the considered 
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models. 
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